Spatially aware deep learning reveals tumor heterogeneity patterns that encode distinct kidney cancer states.

Clear cell renal cell carcinoma (ccRCC) is molecularly heterogeneous, immune infiltrated, and selectively sensitive to immune checkpoint inhibition (ICI). However, the joint tumor-immune states that mediate ICI response remain elusive. We develop spatially aware deep-learning models of tumor and immune features to learn representations of ccRCC tumors using diagnostic whole-slide images (WSIs) in untreated and treated contexts (n = 1,102 patients). We identify patterns of grade heterogeneity in WSIs not achievable through human pathologist analysis, and these graph-based "microheterogeneity" structures associate with PBRM1 loss of function and with patient outcomes. Joint analysis of tumor phenotypes and immune infiltration identifies a subpopulation of highly infiltrated, microheterogeneous tumors responsive to ICI. In paired multiplex immunofluorescence images of ccRCC, microheterogeneity associates with greater PD1 activation in CD8+ lymphocytes and increased tumor-immune interactions. Our work reveals spatially interacting tumor-immune structures underlying ccRCC biology that may also inform selective response to ICI.

Spatially aware deep learning reveals tumor heterogeneity patterns that encode distinct kidney cancer states.

Clear cell renal cell carcinoma (ccRCC) is molecularly heterogeneous, immune infiltrated, and selectively sensitive to immune checkpoint inhibition (ICI). Established histopathology paradigms like nuclear grade have baseline prognostic relevance for ccRCC, although whether existing or novel histologic features encode additional heterogeneous biological and clinical states in ccRCC is uncertain. Here, we developed spatially aware deep learning models of tumor- and immune-related features to learn representations of ccRCC tumors using diagnostic whole-slide images (WSI) in untreated and treated contexts (n = 1102 patients). We discovered patterns of nuclear grade heterogeneity in WSI not achievable through human pathologist analysis, and these graph-based "microheterogeneity" structures associated with PBRM1 loss of function, adverse clinical factors, and selective patient response to ICI. Joint computer vision analysis of tumor phenotypes with inferred tumor infiltrating lymphocyte density identified a further subpopulation of highly infiltrated, microheterogeneous tumors responsive to ICI. In paired multiplex immunofluorescence images of ccRCC, microheterogeneity associated with greater PD1 activation in CD8+ lymphocytes and increased tumor-immune interactions. Thus, our work reveals novel spatially interacting tumor-immune structures underlying ccRCC biology that can also inform selective response to ICI.

Induction of a colitogenic phenotype in Th1-like cells depends on interleukin-23 receptor signaling.

Interleukin-23 receptor plays a critical role in inducing inflammation and autoimmunity. Here, we report that Th1-like cells differentiated in vitro with IL-12 + IL-21 showed similar IL-23R expression to that of pathogenic Th17 cells using eGFP reporter mice. Fate mapping established that these cells did not transition through a Th17 cell state prior to becoming Th1-like cells, and we observed their emergence in vivo in the T cell adoptive transfer colitis model. Using IL-23R-deficient Th1-like cells, we demonstrated that IL-23R was required for the development of a highly colitogenic phenotype. Single-cell RNA sequencing analysis of intestinal T cells identified IL-23R-dependent genes in Th1-like cells that differed from those expressed in Th17 cells. The perturbation of one of these regulators (CD160) in Th1-like cells inhibited the induction of colitis. We thus uncouple IL-23R as a purely Th17 cell-specific factor and implicate IL-23R signaling as a pathogenic driver in Th1-like cells inducing tissue inflammation.

Integrative clinical and molecular characterization of translocation renal cell carcinoma.

Translocation renal cell carcinoma (tRCC) is a poorly characterized subtype of kidney cancer driven by MiT/TFE gene fusions. Here, we define the landmarks of tRCC through an integrative analysis of 152 patients with tRCC identified across genomic, clinical trial, and retrospective cohorts. Most tRCCs harbor few somatic alterations apart from MiT/TFE fusions and homozygous deletions at chromosome 9p21.3 (19.2% of cases). Transcriptionally, tRCCs display a heightened NRF2-driven antioxidant response that is associated with resistance to targeted therapies. Consistently, we find that outcomes for patients with tRCC treated with vascular endothelial growth factor receptor inhibitors (VEGFR-TKIs) are worse than those treated with immune checkpoint inhibitors (ICI). Using multiparametric immunofluorescence, we find that the tumors are infiltrated with CD8+ T cells, though the T cells harbor an exhaustion immunophenotype distinct from that of clear cell RCC. Our findings comprehensively define the clinical and molecular features of tRCC and may inspire new therapeutic hypotheses.

Building Tools for Machine Learning and Artificial Intelligence in Cancer Research: Best Practices and a Case Study with the PathML Toolkit for Computational Pathology.

Imaging datasets in cancer research are growing exponentially in both quantity and information density. These massive datasets may enable derivation of insights for cancer research and clinical care, but only if researchers are equipped with the tools to leverage advanced computational analysis approaches such as machine learning and artificial intelligence. In this work, we highlight three themes to guide development of such computational tools: scalability, standardization, and ease of use. We then apply these principles to develop PathML, a general-purpose research toolkit for computational pathology. We describe the design of the PathML framework and demonstrate applications in diverse use cases. PathML is publicly available at www.pathml.com.

Integrative molecular characterization of sarcomatoid and rhabdoid renal cell carcinoma.

Sarcomatoid and rhabdoid (S/R) renal cell carcinoma (RCC) are highly aggressive tumors with limited molecular and clinical characterization. Emerging evidence suggests immune checkpoint inhibitors (ICI) are particularly effective for these tumors, although the biological basis for this property is largely unknown. Here, we evaluate multiple clinical trial and real-world cohorts of S/R RCC to characterize their molecular features, clinical outcomes, and immunologic characteristics. We find that S/R RCC tumors harbor distinctive molecular features that may account for their aggressive behavior, including BAP1 mutations, CDKN2A deletions, and increased expression of MYC transcriptional programs. We show that these tumors are highly responsive to ICI and that they exhibit an immune-inflamed phenotype characterized by immune activation, increased cytotoxic immune infiltration, upregulation of antigen presentation machinery genes, and PD-L1 expression. Our findings build on prior work and shed light on the molecular drivers of aggressivity and responsiveness to ICI of S/R RCC.

An IL-27-Driven Transcriptional Network Identifies Regulators of IL-10 Expression across T Helper Cell Subsets.

Interleukin-27 (IL-27) is an immunoregulatory cytokine that suppresses inflammation through multiple mechanisms, including induction of IL-10, but the transcriptional network mediating its diverse functions remains unclear. Combining temporal RNA profiling with computational algorithms, we predict 79 transcription factors induced by IL-27 in T cells. We validate 11 known and discover 5 positive (Cebpb, Fosl2, Tbx21, Hlx, and Atf3) and 2 negative (Irf9 and Irf8) Il10 regulators, generating an experimentally refined regulatory network for Il10. We report two central regulators, Prdm1 and Maf, that cooperatively drive the expression of signature genes induced by IL-27 in type 1 regulatory T cells, mediate IL-10 expression in all T helper cells, and determine the regulatory phenotype of colonic Foxp3+ regulatory T cells. Prdm1/Maf double-knockout mice develop spontaneous colitis, phenocopying ll10-deficient mice. Our work provides insights into IL-27-driven transcriptional networks and identifies two shared Il10 regulators that orchestrate immunoregulatory programs across T helper cell subsets.

Checkpoint Blockade Immunotherapy Induces Dynamic Changes in PD-1-CD8+ Tumor-Infiltrating T Cells.

An improved understanding of the anti-tumor CD8+ T cell response after checkpoint blockade would enable more informed and effective therapeutic strategies. Here we examined the dynamics of the effector response of CD8+ tumor-infiltrating lymphocytes (TILs) after checkpoint blockade therapy. Bulk and single-cell RNA profiles of CD8+ TILs after combined Tim-3+PD-1 blockade in preclinical models revealed significant changes in the transcriptional profile of PD-1- TILs. These cells could be divided into subsets bearing characterstics of naive-, effector-, and memory-precursor-like cells. Effector- and memory-precursor-like TILs contained tumor-antigen-specific cells, exhibited proliferative and effector capacity, and expanded in response to different checkpoint blockade therapies across different tumor models. The memory-precursor-like subset shared features with CD8+ T cells associated with response to checkpoint blockade in patients and was compromised in the absence of Tcf7. Expression of Tcf7/Tcf1 was requisite for the efficacy of diverse immunotherapies, highlighting the importance of this transcriptional regulator in the development of effective CD8+ T cell responses upon immunotherapy.

Genome-wide discovery of epistatic loci affecting antibiotic resistance in Neisseria gonorrhoeae using evolutionary couplings.

Genome analysis should allow the discovery of interdependent loci that together cause antibiotic resistance. In practice, however, the vast number of possible epistatic interactions erodes statistical power. Here, we extend an approach that has been successfully used to identify epistatic residues in proteins to infer genomic loci that are strongly coupled. This approach reduces the number of tests required for an epistatic genome-wide association study of antibiotic resistance and increases the likelihood of identifying causal epistasis. We discovered 38 loci and 240 epistatic pairs that influence the minimum inhibitory concentrations of 5 different antibiotics in 1,102 isolates of Neisseria gonorrhoeae that were confirmed in a second dataset of 495 isolates. Many known resistance-affecting loci were recovered; however, the majority of associations occurred in unreported genes, such as murE. About half of the discovered epistasis involved at least one locus previously associated with antibiotic resistance, including interactions between gyrA and parC. Still, many combinations involved unreported loci and genes. While most variation in minimum inhibitory concentrations could be explained by identified loci, epistasis substantially increased explained phenotypic variance. Our work provides a systematic identification of epistasis affecting antibiotic resistance in N. gonorrhoeae and a generalizable approach for epistatic genome-wide association studies.

Induction and transcriptional regulation of the co-inhibitory gene module in T cells.

The expression of co-inhibitory receptors, such as CTLA-4 and PD-1, on effector T cells is a key mechanism for ensuring immune homeostasis. Dysregulated expression of co-inhibitory receptors on CD4+ T cells promotes autoimmunity, whereas sustained overexpression on CD8+ T cells promotes T cell dysfunction or exhaustion, leading to impaired ability to clear chronic viral infections and diseases such as cancer1,2. Here, using RNA and protein expression profiling at single-cell resolution in mouse cells, we identify a module of co-inhibitory receptors that includes not only several known co-inhibitory receptors (PD-1, TIM-3, LAG-3 and TIGIT) but also many new surface receptors. We functionally validated two new co-inhibitory receptors, activated protein C receptor (PROCR) and podoplanin (PDPN). The module of co-inhibitory receptors is co-expressed in both CD4+ and CD8+ T cells and is part of a larger co-inhibitory gene program that is shared by non-responsive T cells in several physiological contexts and is driven by the immunoregulatory cytokine IL-27. Computational analysis identified the transcription factors PRDM1 and c-MAF as cooperative regulators of the co-inhibitory module, and this was validated experimentally. This molecular circuit underlies the co-expression of co-inhibitory receptors in T cells and identifies regulators of T cell function with the potential to control autoimmunity and tumour immunity.

3-Aminobenzamide Blocks MAMP-Induced Callose Deposition Independently of Its Poly(ADPribosyl)ation Inhibiting Activity.

Cell wall reinforcement with callose is a frequent plant response to infection. Poly(ADP-ribosyl)ation is a protein post-translational modification mediated by poly(ADP-ribose) polymerases (PARPs). Poly(ADP-ribosyl)ation has well-known roles in DNA damage repair and has more recently been shown to contribute to plant immune responses. 3-aminobenzamide (3AB) is an established PARP inhibitor and it blocks the callose deposition elicited by flg22 or elf18, two microbe-associated molecular patterns (MAMPs). However, we report that an Arabidopsis parp1parp2parp3 triple mutant does not exhibit loss of flg22-induced callose deposition. Additionally, the more specific PARP inhibitors PJ-34 and INH2BP inhibit PARP activity in Arabidopsis but do not block MAMP-induced callose deposition. These data demonstrate off-target activity of 3AB and indicate that 3AB inhibits callose deposition through a mechanism other than poly(ADP-ribosyl)ation. POWDERY MILDEW RESISTANT 4 (PMR4) is the callose synthase responsible for the majority of MAMP- and wound-induced callose deposition in Arabidopsis. 3AB does not block wound-induced callose deposition, and 3AB does not reduce the PMR4 mRNA abundance increase in response to flg22. Levels of PMR4-HA protein increase in response to flg22, and increase even more in flg22 + 3AB despite no callose being produced. The callose synthase inhibitor 2-deoxy-D-glucose does not cause similar impacts on PMR4-HA protein levels. Beyond MAMPs, we find that 3AB also reduces callose deposition induced by powdery mildew (Golovinomyces cichoracearum) and impairs the penetration resistance of a PMR4 overexpression line. 3AB thus reveals pathogenesis-associated pathways that activate callose synthase enzymatic activity distinct from those that elevate PMR4 mRNA and protein abundance.

Erratum: The neuropeptide NMU amplifies ILC2-driven allergic lung inflammation.

This corrects the article DOI: 10.1038/nature24029.

The neuropeptide NMU amplifies ILC2-driven allergic lung inflammation.

Type 2 innate lymphoid cells (ILC2s) both contribute to mucosal homeostasis and initiate pathologic inflammation in allergic asthma. However, the signals that direct ILC2s to promote homeostasis versus inflammation are unclear. To identify such molecular cues, we profiled mouse lung-resident ILCs using single-cell RNA sequencing at steady state and after in vivo stimulation with the alarmin cytokines IL-25 and IL-33. ILC2s were transcriptionally heterogeneous after activation, with subpopulations distinguished by expression of proliferative, homeostatic and effector genes. The neuropeptide receptor Nmur1 was preferentially expressed by ILC2s at steady state and after IL-25 stimulation. Neuromedin U (NMU), the ligand of NMUR1, activated ILC2s in vitro, and in vivo co-administration of NMU with IL-25 strongly amplified allergic inflammation. Loss of NMU-NMUR1 signalling reduced ILC2 frequency and effector function, and altered transcriptional programs following allergen challenge in vivo. Thus, NMUR1 signalling promotes inflammatory ILC2 responses, highlighting the importance of neuro-immune crosstalk in allergic inflammation at mucosal surfaces.

Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma.

Although human tumours are shaped by the genetic evolution of cancer cells, evidence also suggests that they display hierarchies related to developmental pathways and epigenetic programs in which cancer stem cells (CSCs) can drive tumour growth and give rise to differentiated progeny. Yet, unbiased evidence for CSCs in solid human malignancies remains elusive. Here we profile 4,347 single cells from six IDH1 or IDH2 mutant human oligodendrogliomas by RNA sequencing (RNA-seq) and reconstruct their developmental programs from genome-wide expression signatures. We infer that most cancer cells are differentiated along two specialized glial programs, whereas a rare subpopulation of cells is undifferentiated and associated with a neural stem cell expression program. Cells with expression signatures for proliferation are highly enriched in this rare subpopulation, consistent with a model in which CSCs are primarily responsible for fuelling the growth of oligodendroglioma in humans. Analysis of copy number variation (CNV) shows that distinct CNV sub-clones within tumours display similar cellular hierarchies, suggesting that the architecture of oligodendroglioma is primarily dictated by developmental programs. Subclonal point mutation analysis supports a similar model, although a full phylogenetic tree would be required to definitively determine the effect of genetic evolution on the inferred hierarchies. Our single-cell analyses provide insight into the cellular architecture of oligodendrogliomas at single-cell resolution and support the cancer stem cell model, with substantial implications for disease management.

A Distinct Gene Module for Dysfunction Uncoupled from Activation in Tumor-Infiltrating T Cells.

Reversing the dysfunctional T cell state that arises in cancer and chronic viral infections is the focus of therapeutic interventions; however, current therapies are effective in only some patients and some tumor types. To gain a deeper molecular understanding of the dysfunctional T cell state, we analyzed population and single-cell RNA profiles of CD8(+) tumor-infiltrating lymphocytes (TILs) and used genetic perturbations to identify a distinct gene module for T cell dysfunction that can be uncoupled from T cell activation. This distinct dysfunction module is downstream of intracellular metallothioneins that regulate zinc metabolism and can be identified at single-cell resolution. We further identify Gata-3, a zinc-finger transcription factor in the dysfunctional module, as a regulator of dysfunction, and we use CRISPR-Cas9 genome editing to show that it drives a dysfunctional phenotype in CD8(+) TILs. Our results open novel avenues for targeting dysfunctional T cell states while leaving activation programs intact.